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Abstract 


A numerically efficient algorithm for the fast evaluation of time-dependent potentials/fields, radiated by 
highly elongated three-dimensional (3D) source/observer constellations is presented. This 
computational scheme complements the family of Nonuniform Grid Time Domain (NGTD) based 
algorithms (A. Boag, V. Lomakin, and E. Michielssen, IEEE Trans. Antennas and Propagation, vol. 54, 
no. 7, pp. 1943-1951, July 2006) and (J. Meng, A. Boag, V. Lomakin, and E. Michielssen, J. Comp. 
Physics, vol. 229, no. 22, pp. 8430-8444, 1 Nov. 2010.). At the core of the generic NGTD-approach is 
the idea of using retarded time coordinate and amplitude compensation in order to transform the fields 
radiated by spatially confined and temporally bandlimited sources into slowly varying functions 
amenable to spatial sampling and subsequent interpolation. The proposed algorithm uses the NGTD- 
approach in a multilevel fashion while seeking improved performance, in particular, for elongated 
geometries, by analogy to the time harmonic algorithm described in (N. Costa and A. Boag, 
International Journal of Numerical Modelling: Electronic Networks, Devices and Fields, vol. 25, pp. 
645-655, 2012). To that end, sampling rates are derived for the construction of nonuniform Cartesian 
grids used for representation of the left- and right-propagating fields. Elongated geometries allow for the 
design of exceedingly sparse nonuniform sampling grids comprising roughly a constant (i.e. source 
domain size independent) number of points. The proposed multilevel NGTD algorithm employs such 
compressed sampling representations for both outgoing and incoming fields thus symmetrizing the 
functionality of the sampling grids. This allows us to reduce the per-time-step computational complexity 
(CC) from O(NlogN) of the conventional ML-NGTD (N. Costa and A. Boag, IEEE AP-S International 
Symposium, 2210, Memphis, TN, USA, July 2014) to O(N), with N being the number of spatial source/ 
observer points. Numerical implementation, accuracy, and linear complexity of the proposed algorithm 
will be demonstrated by numerical examples. 






Biography 





ur N 


Amir Boag received the B.Sc. degree in electrical engineering and the B.A. degree in 
physics in 1983, both Summa Cum Laude, the M.Sc. degree in electrical engineering in 
1985, and the Ph.D. degree in electrical engineering in 1991, all from Technion - Israel 
Institute of Technology, Haifa, Israel. 

From 1991 to 1992 he was on the Faculty of the Department of Electrical Engineering at 
the Technion. From 1992 to 1994 he has been a Visiting Assistant Professor with the 
Electromagnetic Communication Laboratory of the Department of Electrical and Computer 
Engineering at the University of Illinois at Urbana-Champaign. In 1994, he joined Israel 
Aircraft Industries as a research engineer and became a manager of the Electromagnetics 
Department in 1997. Since 1999, he is with the Physical Electronics Department of the 
School of Electrical Engineering at Tel Aviv University, where he is currently a Professor. 
Dr. Boag's interests are in computational electromagnetics, wave scattering, imaging, and 
design of antennas and optical devices. He has published over 110 journal articles and 
presented more than 250 conference papers on electromagnetics and acoustics. Prof. Boag 
is an Associate Editor for IEEE Transactions on Antennas and Propagation. He is a Fellow 
of the Electromagnetics Academy. In 2008, Amir Boag was named a Fellow of the IEEE 
for his contributions to integral equation based analysis, design, and imaging techniques. 


W Overview 


Q Motivation - TDDFT 


O Non-Uniform Grid (NG) 
Representation on Elongated 
Domains 


QO Multilevel Algorithm 
a Fast Time Domain Algorithms 
a Conclusions 


% MOTIVATION 


> Time Domain Scattering From Elongated Objects 


> Time Domain Density Functional Theory (TDDFT) 
analysis of elongated 


> Attaining Favorable Computational Complexity 





Scattering by a C,,H, Carbyne chain 


TDDFT — An Introduction 
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wy Maxwell's Equations within TDDFT 


Our time-dependent Kohn-Sham (TDKS) equation of interest is then 
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where ¢, is the usual Kohn-Sham potential 
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The induced vector potential A;,,; depends on the paramagnetic current 
density 





Computation of Potentials 


Gauge-free relations can be expressed as 
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Two gauge-fixing conditions are often used: 
Coulomb gauge: V-A=0 Lorenz gauge: ô + <V -A=0 
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the Coulomb gauge is typically used in TDDFT without the transverse 
vector potential A7. We prefer Lorenz Gauge. 


wW Marching on in Time - TDDFT Solver 


At each time step: 
> Time propagation for N orbitals - y,; (1, t) 


> Evaluation of Source Densities - po(r,t), j(r,t) 
> Computation of Potentials - d&(r,t), A(r,t) 





Computation of Retarded Potentials 


The 4-vector of retarded potentials is evaluated: 
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using 4-vector of sources and smoothed discrete kernel: 
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The sources at retarded times are interpolated: 


S(r,.r)= %, wi'-t,)S (r,,t,) 
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Per time step Computational Complexity (CC) of Direct Evaluation: O(N?) 


Sampling of Far Zone Potentials 


Consider the far zone potential of sources confined to a subdomain V 
centered at r = (0,0, Z) 
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We define a delay and amplitude compensated potential: 
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Right propagating wave 






with retarded time T =t—F/c 
and corrected distance 
to the center r= Jlr—rl 42h? 


T =t—z/c= const. 
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Sampling of Far Zone Potentials 


Delay and Amplitude Compensation: 
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Spatial Interpolation: 
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Delay and Amplitude Restoration: 
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Sources are bandlimited to ,,,, 


Sampling criteria can be derived in frequency domain 


Frequency Domain Potentials 


Consider the far zone potential of sources confined to a subdomain V 
centered at r 
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We define a phase and amplitude compensated potential: 
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with the phase correction kr 


and corrected distance to the center 7 = ./Ir—r I? +2h? 


Spatial Sampling / Non-uniform Grid 


Consider the phase of the compensated potential: 
O(r) =kF-Ir—r'l) = k( vir FP +h? -Ir—r' ) = 
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We define local frequencies with respect to the Cartesian coordinates as: 
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Then sampling rates are estimated based on local Nyquist criterion: 
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Spatial Sampling / Non-uniform Grid (NG) 


The sampling rates: 
ff, =Q,(2kh* /74+C,), aeFZV) & Q,>1 
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The sampling rates for bandlimited TD signal (K,,,, = @nax /©): 
f, =Q,(2k,,,.0 /2#+C,), a@EFZV) & Q,>1 


Fi, (@)=Q, ,(Qg+ Dlalky hl +C,,), 2EFZV) & Q,,>1 


Number of sampling is bounded — O(1)! 


ML-NGTD- Preliminary Stage 


Hierarchical Multilevel Domain Decomposition: 
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Finest level subdomains 
At finest level of division: R, S Ain 


We assume identical half-width A for all subdomains. 


% ML-NGTD - Near Field Contribution 


At level L 
Yyro= > Worn rev, 
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Near Field Computational work — O(N) 
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ML-NGTD — Far Field Initialization 


At level L 
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Far Field Initialization Computational work — O(N) 


ML-NGTD — Outgoing Field Aggregation 


Transition from level / to /-1 
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Repeated from /=L up to /=2 
Computational work — O(N+ N/2+N/4+...) = OW) 


ML-NGTD - Incoming Field Initialization 


At level /=2 
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Incoming Field Initialization Computational work — O(1 ) 


ML-NGTD - Incoming Field Aggregation 


Transition from level /-1 to / 
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Repeated from /=3 up to /=L 
Computational work — O(...+N/4+N/2+N) = O(N) 


ML-NGTD - Final Field Aggregation 


At level L 
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Computational work — O(N) 
Total Complexity — O(N)! 
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wy Complexity of TD Algorithms 


> FFT-based methods — O(N[logN +log?N,]) 

> Plane wave time domain (PWTD) — O(Nlog?N) 

> Non-uniform grid time domain (NGTD) — OWVlogN) 

> Cartesian non-uniform grid time domain (NGTD) — O(Nlog? N) 


> NGTD for elongated geometries — O(N) 


% Conclusions 


o A novel Multilevel Cartesian non-uniform grid time domain (NGTD) 
algorithm that permits fast potential/field evaluation with Linear CC 
for given source constellations in elongated geometries was 
presented. 


a This scheme demonstrates the geometry adaptivity of the generic 
NG approach. 


a Implementation and integration with TDDFT should be useful in 
accelerating first principles solution for elongated nano-antennas 


a Also MOT solvers for conventional linear arrays of antennas and 
particle chains can be accelerated. 
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